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1. Introduction 



In a recent paper M one of the present authors has proposed a numerical implementation 
of a method of successive approximations that allows automatic construction of the 
order reduction that contains the physical, non-runaway, solutions of a large class of 
singular differential equations, including the classical equation of motion of charged 
particles with radiation reaction |2| and fourth-order equations appearing in theories 
of gravitation with a quadratic lagrangian || and in the study of quantum corrections 
to Einstein equations @|. Apart from its practical interest, the convergence of the 
numerical method provides indirect but convincing evidence of the convergence of the 
analytic method. 

The goal of this letter is twofold: we want to discuss the phenomenon of 
preacceleration in the frame of the method of successive approximations and to produce 
a more physical exact example in which the method can be analyzed in full detail. 



2. Preacceleration 



Among the puzzling properties of the Lorentz-Dirac equation |J, which describes the 
motion of a radiating point charge, the preacceleration is one of the consequences of its 
singular structure. Let us consider the non-relativistic approximation to the Lorentz- 
Dirac equation (the so-called Abraham-Lorentz equation) in the case of a charge e that 
moves in a straight line under the action of an external force per unit mass /(£): 

X = f(t)+T X, (1) 

where 

2e 2 

T ° S <W- (2) 
It is well known that if one eliminates the runaway solutions, the physical motion is 
described by the integro-differential equation |5|, |2], [L| 

x = \ e- u f{t + T u)du. (3) 
Jo 

If the external force is f(t) = foS(t) it is enough to insert this expression in 
equation (f|) to obtain the order reduction 

x{t) = I r (4) 
{ if t > 0. 

According to this the charge would start accelerating before the pulse reaches it. 

This phenomenon has been widely discussed in connection with the smallness of To 
and taking into account the limitations of the classical theory, but we want to analyze 
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it here from the point of view of the method of successive approximations [|TJ, which 
starts from the approximation that neglects completely the radiation reaction 

x = e (t) = f(t), (5) 

and iteratively constructs approximate reductions by substituting the previous 
approximation on the right hand side of P) 

n+l 

x = e n+1 (t) = f(t) + r e' n (t) = r k f {k \t). (6) 

fc=0 

Under the appropriate mathematical conditions, this method will converge to the exact 
reduction 

oo 

x = Q(t) = ^^f {k \t), (7) 

fc=0 

which is precisely the Taylor expansion of (|3|). 

Now, one of the main hypothesis in the method of successive approximations is 
Bhabha's remark || that the physical solutions are precisely those that are regular in 
the limit r — > 0, where according to (|3]) and (0) one recovers the second order equation 



x = f(t). (8) 



But we can see that the preaccelerated solution ([|) is divergent in the limit r — > 0, 
exactly as the remaining pathological solutions. In consequence, it cannot be constructed 
by the (analytical or numerical) method of successive approximations. One might see 
this as a limitation of the latter method, but we think that it is rather a limitation of 
the Abraham-Lorentz and similar equations. 

To make clear our point of view we will consider a pulse of small but non null width. 
For simplicity we will take a Gaussian pulse, 

/( t) = -4=e-^ 2 , (9) 

but one could also consider any other pulse that recovers in the limit e — > the value 
fo5(t). After inserting (||) in (|3]) one gets the Newtonian equation of motion that 
contains the physical solutions: 

■ ± = A e */ro e eV^ erfc (t + J_Y (10) 
2r \e 2tJ k j 

This would be precisely the reduction constructed by the method of successive 
approximations. By using the properties of the complementary error function, it is 
easy to see that in the limit r — > one recovers the radiationless result (^]), while for 
e — > one gets the preaccelerated solution We see, thus, that these two limits do 
not commute. 
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In our opinion the delta function obtained in the limit e —>■ is beyond the field 
of applicability of Lorentz-Dirac and Abraham-Lorentz equations, for which one has to 
assume that the applied force and acceleration do not change too much across a time 
interval of length tq, i.e. the radiation reaction cannot be too important along such a 
tiny interval. Since this assumption is not met by the delta function, this often useful 
limit is not applicable here and one has necessarily to consider pulses of width larger 
than tq. This opinion is in agreement with the point of view of authors || |7|, Q that 
stress that the analyticity with respect to r is a fundamental hypothesis, which is used 
in standard derivations of the Lorentz-Dirac equation. 



3. An exact example 



In reference Jl| we discussed a linear one- dimensional exact example in the frame of the 
method of succesive approximations. We want to analyze now a three-dimensional exact 
example that, though still linear, has a clearer physical meaning and will contribute to 
our confidence on the convergence of the method of successive approximations under 
appropriate conditions. Let us consider a charge e that moves in an external magnetic 
field B, as happens in some astrophysical contexts |5| or in particle accelerators [IU|. In 
the non-relativistic approximation the equation of motion is 



x = Q x x + r x, 
where we have introduced the cyclotron frequency 

m 



(12) 



which we will assume to be uniform and constant. Starting from the lowest order 
approximation 



x = ©o = CI x x, 

we can construct successive approximations by using repeatedly 

~d@ r 



® n+ l = X X + Tq 



Ot 



+ (X • V x ) n + (0 n • Vi) @ n 



It is straightforward to check that the successive approximations are 
n = a n fl x x - /3 n x_L, 

where 



fi X 



(13) 



(14) 



(15) 



(16) 
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is the component of the velocity perpendicular to the magnetic field and the constant 
coefficients are given by the recurrence 

a n+1 = 1 - 2r a n p n , (17) 
f3 n+1 =r (n 2 a 2 n -(3 2 n ), (18) 



and the initial conditions «o = 1 and (3q = 0. 

Recurrence (18) has two fixed points P± = (a±, j3±) with 



^= ± V 2r n ^ 




P± = 5 > (20) 

2r 

but a linear stability analysis proves that P_ is always unstable and that P + is 
asymptotically stable for 



To n < v 3 + 2 ^ w 64 (21) 

Furthermore, a simple bifurcation diagram in the dimensionless variables (a n , ro/3 n ) 
shows that the initial condition (1,0) is in the basin of attraction of P + and, in 
consequence, that the method of successive approximations will in fact converge in 
the range (pi] ) to the Newtonian equation 

x = a + Q, x x — (3 + xj_, (22) 

which contains precisely the physical (non runaway) solutions for x found by Plass ||. 
Notice that all the approximations (|I5|), as well as the exact order reduction (|22|), 
are orthogonal to the magnetic field, and that the exact reduction exists even when 
the method fails. This is not surprising because most approximation methods have 
limited ranges of applicability Moreover, in this case the method will converge in all 
practical situations because the cyclotron frequency is always very small compared to 
1/tq. This simple but illustrative exact example reinforces our conviction that the 
numerical approximation method |l| will converge in many cases of interest. 
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